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Abstract 

We investigate the possibility to use atom interferometers to detect 
gravitational waves. We discuss the interaction of gravitational waves 
with an atom interferometer and analyze possible schemes. 
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1 Introduction 

The direct detection of Gravitational Waves (GWs) is one of the most exciting 
scientific goals because it would improve our understanding of laws governing 
the universe and provide new means to observe it. The most sensitive detectors 
which are already operating, under construction or being planned are based on 
optical interferometers [Tl[21[3]. In most cases, however, the sensitivities are 
only marginally sufficient to detect the expected signals, detectors have large 
sizes ranging from a few km on the Earth (Virgo, LIGO) to millions of km in 
space (LISA) and the operating frequency ranges are limited. Therefore it is 
of great interest to investigate alternative schemes that can lead to a higher 
sensitivity, smaller sizes and extend the frequency range of the detectors. 

In recent years, matter-wave interferometry with neutral atoms has under- 
gone an impressive development due to the increasing ability to control the in- 
ternal and external atomic degrees of freedom using laser manipulation methods 
[H [51 m 13 • Atom Interferometers ( AIs) are already competing with state-of-art 
optical interferometers in terms of sensitivity. This was demonstrated experi- 
mentally for gravity acceleration 0, gravity gradients [S], inertial and rotation 
effects [inilll]. Other experiments, planned or presently in progress, to inves- 
tigate properties of gravitational field by AI concern accurate measurements of 
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G [m [IH] , test of the equivalence principle [2] , detection of the Lense-Thirring 
effect fTS], deviations from the 1/r^ Newtonian law for small distances [16] . 

In analogy to optical interferometers, in atom interferometers atomic wave 
packets are split and recombined giving rise to an interference signal. Different 
schemes can be used for splitting, reflecting and recombining the atoms. In a 
particular class of interferometers, which is the one relevant in this paper, the 
separation of the atoms is achieved by inducing a transition between internal 
states of the atoms by an electromagnetic field. The spatial separation in this 
case is induced by the momentum recoil and the internal and external states 
of the atoms become entangled. Another approach is to use material gratings. 
This raises however different problems, both conceptual and technical, such as 
the realization and handling of the required nano-structures and it will not be 
considered here. 

In this paper, we discuss the possibility to use AI to detect gravitational 
waves. The interaction between matter waves and gravitational waves was al- 
ready investigated in [TTl [TSl (THl [20] • Recently, due to the experimental ad- 
vances, the interest was revived [HI [HI [23] [211 [15]. The aim of the present 
paper is to analyze possible schemes for interferometers using light fields as 
atom optics components. Compared to |25| . where only the Linet-Tourrenc 
contribution is considered in the eikonal approximation, we take into account 
all the contributions in phase difference inside the interferometer. 

The paper is organized as follows. In Sect. [H we recall the ABCD formalism 
for matter waves and apply it to the calculation of phase shift in atom interfer- 
ometers in the specific case of weak gravitational field when the Hamiltonian is 
at most quadratic in coordinates and conjugate momenta. A detailed discussion 
about Einstein and Fermi coordinates is presented in appendix. We apply the 
results in Sect. [3] to derive the phase shift signal for two specific atom inter- 
ferometer configurations. Finally, in Sect. [4] we discuss possible experimental 
schemes and evaluate the sensitivity for the detection of gravitational waves. 

2 The ABCD matrices for matter waves and 
phase-shift formula for atom interferometers 

In this section we recall the ABCD formalism for matter waves and apply it to 
calculation of phase shift formula. The discussion is based on the relativistic 
Schroedinger-type equation for atom waves and its analysis in |26| and references 
therein. 

2.1 The ABCD matrices for matter waves 

In the following, we assume that the Hamiltonian relative to the motion of the 
center of mass is a quadratic polynominal of momentum and position operators. 
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as in most cases of relevance in AI: 
1 
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where a, j3, 7, 6 are suitable square matrices coefRcients for the quadratic terms 
(with 5 = —a, where the tilde indicates the transposed matrix); g is the grav- 
ity vector field and / is an external vector field. M* is the relativistic mass 
(M* = Mq/a/I — w^/c^, where Mq is the rest mass). The evolution of the 
wave packets by this Hamiltonian, via the Ehrenfest theorem, can be obtained 
through Hamilton's equations 
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The integral of Hamilton's equations can be written through the ABCD matrices 
as 
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with T the time ordering operator. A perturbative expansion leads to [27j : 

Mit,to) = l+ f Tit')dt' + f dt' f rit')Tit")dt" + ... (8) 
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The Eq. ([7]) can be used to find the ABCD matrices which determine the 
evolution of the wave packets in presence of the GW, h. If we consider the 
simple case in which the GW is the only (weak) field, in the two coordinate 
systems discussed in Appendix A, for a single Fourier component we have up 
to the first order in h{^l)■. 



Fermi coordinates 
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where V is the Minkowskian matrix. 
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2.2 The phase shift formula for atom interferometers 



The total phase difference between the two arms, s and i, of an atom inter- 
ferometer can be expressed as the sum of three terms: the difference in the 
action integral along each path; the difference in the phases imprinted by the 
beam-splitters on the atom waves; the contribution from the splitting of the 
wave packets at the exit of the interferometer [26] : 

1 ^ 

N 

+^\PsD{q- Qsd) - PiD{q~ qio)] (13) 

where Ssj — Ss{tj^i,tj) and Sij — Si{tj^i,tj) are the action integrals along 
s {i) path; ksj{kij) is the momentum transferred to the atoms by the j-th 
beam-splitter along the s (i) arm; q^j and qij are the coordinates of the beam 
splitter/atom interaction; ujsj{ujij) is the angular frequency of the laser beam; 
'dsji'&ij) is the phase of the laser beam at the j-th interaction with the atom; D 
is the exit port. 

Assuming the same input point for the two arms and using the " mid point" 
property in integrating over the space at the output, the phase shift difference 
A(/3 between the two arms (s,i) for an interferometer with N beam splitters can 
be written as: 

Af , Af 



Af Af 



+ (14) 



where Msj{Mij) is the mass of the atom in the s (i) arm and Tj is a proper time 
at the j-th interaction. 



3 Phase calculation for different AI geometries 

In this section, we use the ABCD formalism to find the resulting phase shift A(p 
for typical atom interferometer schemes [1]. The approach is in the frequency 
space of complex Fourier transform in order to describe both the amplitude and 
the phase of the resulting Aip. We consider here the trapezoidal interferometer, 
first suggested in [28], and the parallelogram-shaped interferometer, both in 
Fermi and Einstein coordinates, retaining only terms up to the first order in h. 
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3.1 Trapezoidal AI 

The scheme of a trapezoidal interferometer is shown in Fig. [TJ The interaction of 
the atom with two counter-propagating pairs of copropagating beams (BS1-BS4) 
gives rise to a trapezoidal-shaped closed circuit. By using Eq. (|14p . we obtain 
the expression of the phase shift for the two arms (s,i) of the interferometer: 

AV5 = kiqi + ^ks2{qs2,b + qi2.a) + ^ks3{qt3,a + qs3,a) + ^hiQsAM + qii.a) 



(15) 



where 51^ is the laser frequency and is the frequency of the atomic transition 
involving ground (a) and excited (b) states; di are the proper laser phases. Ex- 
pressing all the qj coordinates through the ABCD matrices, the phase difference 
at the output of the interferometer is then given by: 
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where and are the masses of the atom in the ground and excited state, re- 
spectively; the expression for the A, B, C, D matrices depends on the coordinate 
system. 

Let's first consider Fermi coordinates. From Eq. and Eq. pUj) . wc obtain 
for the phase difference for this configuration 



1 
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Considering Einstein coordinates, from Eq. and Eq. (H^), we obtain 
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where 'diE are the laser phases in the Einstein coordinates system. They are 
different because in Einstein coordinates the index of refraction for the vacuum 
is varying with h(t) or, which is the same, we have an extra (Fourier trans- 
formed) contribution Sk ^ k[h{t, fl) /2]eKp(iQt) to the momentum transfered at 
the beam-splitter positions as a consequence of the apparent photon velocity 
V ^ c[l + {h/2)]. This can be accounted for in the phase terms leading to an 
extra term —Sk{t,fl) in Eq. ([15]), as in spatial beam-splitters '26^. By inserting 
these laser phases and using the coordinates transformation rules in GR between 
the two systems here considered (see Appendix B), the resulting phase shift dif- 
ference is coincident with the one obtained in Fermi coordinates, as expected 
for a scalar quantity which is a physical result in spite of different descriptions. 



3.2 Parallelogram-shaped AI 

The scheme of a parallelogram-shaped AI is shown in Fig. [2l In this case, the 
interaction of the atom with four copropagating laser beams gives rise to a 
parallelogram-shaped closed circuit. 

The phase difference at the output of this interferometer is given by: 

= fci<7i [1 - 2A{T, 0) + A{2T, 0)] + 

4[B(2r.0,-2B(r,0)l(^ + ^ + 5|)+ (19) 

hk^ 

+— f B(2T, T) [D{T, 0)-l]e + §,-^2-^3+^4 
1Mb 

up to the first order in e, where e = [Mb — Ma) /Ma- The only difference from 
the case of a trapezoidal interferometer (Eq. (|16p ) is in the recoil term, which is 
proportional to the relative energy difference between ground and excited states. 
Using Fermi coordinates, from Eq. ([TO]) we obtain: 

' 2 ^\Ma Mb Mb 
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Considering Einstein coordinates, we obtain: 
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The same considerations of paragraph [3A] apply in this case about the iden- 
tity of the results in the two descriptions. 



4 Possible schemes and expected sensitivities 

The results in previous section provide the phase shift at the output of the atom 
interferometer induced by a gravitational wave with amplitude h and frequency 
Q for the typical schemes considered. In order to determine the sensitivity of 
the interferometer, that is the minimum detectable amplitude h{Q), we assume 
shot-noise-limited detection of the atoms, corresponding to a phase noise given 

by Aip{D,)) — rj/V^ where rj is a detection efficiency and TV is the atoms flow 
at the detector. The resulting sensitivity h{^l) (at S/N ratio equal 1) can be 
written as 

where E is a scale factor and f{fl) is the resonance function of the interferometer. 
Neglecting the clock and recoil terms, from Eq. (fT7)) we obtain: 

h{n) = -1- (23) 
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where L is the characteristic hnear dimension of the interferometer, px = mvx 
with vt the transversal velocity acquired by the atoms in the splitting process 
and 2/m = 1/Ma + l/Mt, and where 



From Eq. (|23p . it is evident that in order to achieve the required sensitivity 
while keeping a sufficiently large detection bandwidth it is necessary to realize 
large values of L and pr ■ 

In order to evaluate the performance of these new detectors for GWs, we 
analyzed a few specific cases. It is important to notice that in this analysis we 
did not treat other noise sources that, as in optical GW detectors, can affect 
the performance of the AI detector. Examples are the suspension of optics 
required for the manipulation of the laser beams or the phase noise of the laser 
itself. Based on the work for optical GW detectors and progress in ultrastable 
lasers for future optical clocks, suitable laser sources and suspension systems 
can be envisaged. A detailed analysis of the overall noise budget, including 
technological aspects, is beyond the scope of the present paper. 

Let us consider first an atom interferometer based on a fast beam of hydrogen 
atoms. If we take T — 10^^ s and a length L = 10^ m, similar to present optical 
interferometers detectors, we have vl = 10^ m/s. As shown in Fig[31 in this 
case a sensitivity h{fl) of about 10~'^^ /Hz^^"^ is achieved for vt ~ 10 m/s with a 
flux of lO""^^ atoms/s in the atomic beam [29^. The recoil velocity for a hydrogen 
atom absorbing a Ly — a photon is Vrec = 3.3 m/s. Although the absorption of 
a photon followed by spontaneous emission destroys coherence and cannot be 
used to deflect atomic trajectories in an interferometer, it is conceivable to use 
two-photon Raman transitions between the two hyperfine levels of H ground 
state. A single Raman pulse transfers a velocity vt = 2vrec- Raman transitions 
have already been used in AI based on alkali atoms [1] and the possibility to use 
multiple Raman pulses sequences to increase the enclosed area and the resulting 
sensitivity was also demonstrated [30j . A practical limitation at present would 
be the required power (w 10 W) of laser radiation at the Lya wavelength. 
This is orders of magnitude larger than what can be presently achieved as cw 
radiation |31j but closer to what is produced in pulsed mode 02] ■ An alternative 
scheme is the excitation of a two-photon transition from the ground state to 
long lived excited states [551 [51] . A large recoil can be transfered by combining 
Is — 2s excitation with optical transitions from the 2s state to high lying p 
states [5S]. Such a scheme is also compatible with a ground based apparatus 
because of the negligible vertical displacement of the atomic beam during the 




(24) 
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short total time of flight. Other cases we considered for fast beams of heavier 
atoms do not meet the requirements for the scheme we considered because of the 
difficulty to transfer a large enough transverse momentum to the atoms while 
keeping T small enough in order to keep a large bandwidth. An improvement 
in sensitivity, at the expense of a reduced bandwidth, could be achieved by 
increasing T and correspondingly the linear dimensions of the interferometer. 
In this case, however, a gravity-free apparatus in space should be considered 
(FigEl). 

A different case we considered is an interferometer based on cold atoms. In 
this case, it is more useful to rewrite Eq. in the form: 

(25) 

^ \PtJ niL^ f{n) 

where m is the atomic mass. It is apparent that in this scheme, by relaxing the 
constraint on T, the sensitivity is better the larger is the value of m and the 
smaller is pL- As an example, if = 1 m/s = 2vt , i = 50 m and m « 10^ 
a.m.u., a sensitivity of about 10^^^ /Hz^^^ results at frequencies around 10 mHz 
(Fig. [3]). The long time of flight T w 50 s, for which a gravity-free scheme would 
be required, leads of course to a narrower bandwidth. 



5 Conclusions 

We investigated the possibility of detecting gravitational waves using atom in- 
terferometers based on light fields as beam-splitters. The phase shift at the 
output of the interferometers was calculated for presently known schemes using 
both Einstein and Fermi coordinates. Considering sensitivities of the same or- 
der of magnitude as the ones of present optical gravitational wave detectors, we 
estimated the resulting values for relevant parameters. The results show that 
dedicated technological developments would be needed to achieve the required 
values which are beyond those presently available. New schemes for atom in- 
terferometers, beam splitters, and high flux coherent atomic sources can lead 
to an increase in sensitivity and make atom interferometers competitive with 
other gravitational wave detectors. 



Acknowledgments 

The authors acknowledge useful discussions with Ch. J. Borde. 



10 



Appendix A: Einstein vs Fermi coordinates 



In general relativity (GR) all the coordinates systems are a priori equivalent. 
The predicted physical results do not depend on the specific coordinates sys- 
tem although different descriptions depend on coordinates systems. Generally 
speaking, a change of coordinates is defined by any set of functions 

a;"=a;"(/)^zj"=2/"(a:^) a,/? = 0,1, 2, 3 (26) 

where the invertibility is guaranteed if and only if 

In the following, we refer to the case in which deviations from the Minkowsky 
space of special relativity are due only to GWs in the weak field approximation, 
that is 

gtiu{x) 77^1/ + h^u{x) (28) 

with \h^^\ « 1. We assume that no other field is present. In this case, the 
linearized Einstein field equations admit a plane wave solution for hf^^. It is 
always possible to choose Gaussian synchronous coordinates |36j in which 

50m =.9pO- (1,0,0,0) (29) 

Let's choose the particular coordinates system in which two particles A and B, 
at rest in Minkowsky system, remain at rest even in presence of a GW. These 
are called Einstein coordinates (EC) and the related metrics can be written as 

ds^ = c^dt^ - (1 - h+)dx^ - (1 + h+)dy'^ - dz^ + 2h^dxdy (30) 

for a wave propagating along the z axis; /i+ = h+{ct—z) and hx — hx (ct—z) are 
the amplitudes of the two polarizations states. The world-lines of a free particle 
are geodesies I37| . It is important to note that the proper distance between two 
particles A and B always at rest in this system is varying with the amplitude of 
the GW [SH]: 

d^ = dl + h,, (t) ix], - x^Xx's - x^^) (31) 

where (i§ = (x^ — x\){x'g — x\) in the hypothesis that the particles A and B are 
close enough to consider hij depending only on t. From another point of view, 
we can say that the fiight time of a photon from A to B and back is varying, 
or that we have an index of refraction of the vacuum which is varying with the 
perturbation hij{t). 

The Einstein coordinates are formally the most convenient to describe plane 
GWs; they can be considered as a "wave system (TT gauge)". This is not an 
intuitive system, however, for measurements in a laboratory; as an extension 
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of "classical" approach, indeed, we search for inertial systems in which it is 
possible to preserve the Newtonian idea of "rigid stick" and related measurement 
method. Fermi coordinates (FC) are the best approximation to such a " Galilean 
system" [371 132] ■ We choose a " fiducial observer" (free falling observer) at rest 
in the origin; the system is related to the geodesic line x = y = z = 0. This 
is a Minkowskian system if we disregard small terms as h+ and hx [lOj- The 
spatial axis are built as locally orthogonal coordinate lines whose direction can 
be checked by gyroscopes [321 I2H1 HO]- The transformation between the two 
systems EC and FC are [3511101: 

X = x - ^h+x+ ^h^y 

Y = y+^h+y+^hxx (32) 
Z = z + -!-/i+(x^ - y^) + -^hxxy 

T^t+ -^h+{x^ - y^) + ^hxxy 

where X, Y, Z, T are Fermi coordinates, x, y, z, t are Einstein coordinates, and 
the dot indicates the time derivative. 

Two particles A and £?, initially at rest in FC, move approximately as 

Xa,B = XA,B - -^h+XA.B + -jhy.yA,B 

Ya.B ^ yA,B + ^h+yA^B + ^hxXA,B (33) 
Za,B = ZA,B 

In this case, the proper distance does not change while the distance between the 
two particles does: the index of refraction of the vacuum is 1. Time of fiight of 
a photon between two test masses A and B is the same in both systems: it is a 
physical result, indeed, in spite of the different descriptions. 

It is to be noted that the particle A moves with respect to the particle B as 

subjected to the "tidal" force Fa^i = \mAX^ ■ 

From Eq. ([3^ , it appears that in the study of interaction of GW with exper- 
imental devices, Einstein coordinates are not the most suitable because of the 
very complex motion resulting for an observer in the laboratory frame; Fermi 
coordinates can indeed be considered as the natural extension of a "Cartesian 
inertial system" of the local observer [411 Il2l US] . 



Appendix B: Gauge invariance 

Demonstrating the invariance of the results in Sect. 3 under general gauge trans- 
formations is the subject of ongoing work; here we restrict the discussion to the 
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so called " long wavelength approximation" [J^l . 

It is easy to see that the transformation matrix S from EC to FC of Appendix 
A behaves as 

S = I + 0{h) (34) 

where I is the identity matrix; furthermore from Eql321 we get 

T = t + 0{hL^/c\) (35) 

From Eql34l using Eql35]in the approximation L/A ^ after the insertion of 
the 5k term for proper laser phases, the identity of results in both coordinates 
systems used in Sect. 3 follows. 

It is worth noting that, starting from general FC, it is possible to build a 
simpler " laboratory frame" , that is a rigid coordinates system [23[ I40j , which 
preserves the FC properties in the hypothesis of constant z (wavefront of the 
gravitational plane wave) near the Z — plane (the plane of the interferometer) 
[25., 40,, 44J . Considering for simplicity only the + polarization, the transforma- 
tion law from the EC with metric 

ds^ = c^dt^ - (1 - h+)dx^ - (1 + h+)dy'^ - dz^ (36) 

to the rigid system, is 

t = T 

a; = A + \h+X 

y^Y- ^h+Y (37) 
z^Z 

with transformation matrix 

/ 1 

V 

and metric 

= c^dT^ - dX^ - dY^ - dZ^ + (39) 

+ — {XdXdZ - YdYdZ - cXdXdT + cYdYdT) 
c 

Writing proper ABCD matrices for the rigid coordinates system by using Eq. l39[ 
the same results as in Eq. [TT]or Eq. [20]can be obtained, thus demonstrating the 
identity of the results in rigid, Fermi, and Einstein coordinates systems. 






l-\h 







X 



(38) 
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Detection 



Figure 1: Scheme of the trapezoidal interferometer. Dotted arrows represent 
laser beams acting as beam-splitters (BS1-BS4); bold continuous arrows show 
the relevant momentum transferred to the atom; a: ground internal atomic 
state; b: excited internal atomic state; k: transferred momentum (in h units). 
In the text Ti = T2 = T and T' ^ . 




Detection 



Figure 2: Scheme of a parallelogram interferometer. Dotted arrows represent 
laser beams acting as beam-splitters (BS1-BS4); bold continuous arrows show 
the relevant momentum transferred to the atom; a: ground internal atomic 
state; b: excited internal atomic state; k: transferred momentum (in h units). 
In the text Ti = T2 = T and T' ^ . 



17 



1E-15 

1E-16 

1E-17 

^ 1E-18 
X 

1E-19 

cr 
W 

^ 1E-20 t- 
1E-21 
1E-22 

1E-23 y- 

1E-24 



1E-3 



0.01 



0.1 



10 



QJ2k (Hz) 



100 



1000 



Figure 3: Sensitivity curves of GW detectors based on atom interferometry for 
the three parameter sets discussed in the text: L ~ 10'^ m, vl = 10^ m/s 
(continuous line); L = 2 ■ 10^ m, vl = 5 ■ 10^ m/s (dotted line); L = 50 m, 
Vl = 1 m/s (dashed-dotted line). 
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